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Abstract 

We construct a class of nuclear equations of state based on a schematic potential model, that originates from 
the work of Prakash et. al. [1], which reproduce the results of most microscopic calculations. The equations of 
state are used as input for solving the Tolman-Oppenheimer-Volkov equations for corresponding neutron stars. 
The potential part contribution of the symmetry energy to the total energy is parameterized in a generalized 
form both for low and high values of the baryon density. Special attention is devoted to the construction of the 
symmetry energy in order to reproduce the results of most microscopic calculations of dense nuclear matter. The 
obtained nuclear equations of state are applied for the systematic study of the global properties of a neutron 
star (masses, radii and composition). The calculated masses and radii of the neutron stars are plotted as a 
function of the potential part parameters of the symmetry energy. A linear relation between these parameters, 
the radius and the maximum mass of the neutron star is obtained. In addition, a linear relation between the 
radius and the derivative of the symmetry energy near the saturation density is found. We also address on the 
problem of the existence of correlation between the pressure near the saturation density and the radius. 
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1 Introduction 

Neutron stars (NS) are some of the densest manifestations of massive objects in the universe which provide very 
rich information for testing theories of dense matter physics and also provide a connection among nuclear physics, 
particle physics, statistical physics and astrophysics [ll[31lll[Sl[Sl[71[51[ni [111 EI]- The global aspects of neutron 
stars, such as the masses, radii and composition are determined by solving the so-called Tolman-Oppenheimer- 
Volkov (TOV) equations [121 113j . However there are large variations in predicted radii and maximum masses 
because of the uncertainties in the nuclear equation of state (EOS) near and mainly above the saturation density 

[a[Tl[l3[li[ni[Ti[Ta[20l[2Tl[2l[23[22^ Xhe total energy 

of neutron rich matter (the case of a neutron star) can be written as a sum of two parts. The first one is the 
contribution of the symmetric nuclear matter (which is well known) and the second is the symmetry energy (SE) 
which still is uncertain although several constraints exist from ground state masses (binding energies) and giant 
dipole resonances of laboratory nuclei. A consequence of this uncertainty is that different models predict up to 
a factor of 6, variations in the pressure of neutron star matter near Us, even though the pressure of symmetric 
matter is better known, being nearly zero at the same density. This pressure variation accounts for the nearly 50% 
variation in the predictions of neutron star radii [51 [37] . 

In general, the value of the SE at nuclear saturation density and mainly the density dependence of the SE 
are both difficult to be determined in the laboratory. The motivation of the present work is to propose a new 
parameterization for the potential part of the symmetry energy Egymin) in order to be able to reproduce the 
results of a variety of microscopic models both in low and high values of the baryon density. Especially the trend 
of the symmetry energy just above the equilibrium density ris is a critical factor in determining the neutron star 
radius. 

In order to calculate the global properties of neutron stars (mass, radius ets.) the hydrostatic equilibrium 
equations of Tolman, Oppenheimer and Volkov have to be solved once the equation of state is specified. However, 
the composition of a neutron star still remains uncertain and the construction of the EOS, which is based on the 
ingredients of the NS's and the kind of interactions which characterize them, is subjected to several assumptions. In 
any case the calculated EOS has to satisfy the following requirements [55]: i) It must display the correct saturation 
point for symmetric nuclear matter (SNM); ii) it must give a SE compatible with nuclear phenomenology especially 



at high densities; iii) for SNM the incompressibility at saturation must be compatible with the values extracted 
from phenomenology; iv) both for neutron matter and SNM the speed of sound must not exceed the speed of light 
(causality condition), at least up to the relevant densities. 

In the present work we consider that the neutron star core is composed only by an uncharged mixture of 
neurons, protons and electrons in equilibrium with respect to the weak interaction (/3-stable matter). However in 
general, in the range of densities n > Us, the hadronic phase of superdense matter with a rich spectrum of particles 
(hyperons, baryonic resonances, 7t~ and K~ mesons, and a small portion of leptons) is realized. The model which 
is used for the construction of the EOS is a generalization of a schematic potential model based on a previous work 
of Prakash et al. [T] . The model reproduces the results of most microscopic calculations of dense matter [25] . It is 
worthwhile to notice that there many ways to determine the equation of state through the many-body approach 
of interacting hadrons. Some of the most recent ones are based on variational methods [351 IMl HO] and some are 
based on microscopic calculations |41| . In order to face the problem that stems from the uncertain behavior of 
the SE at high densities we perform a suitable parameterization both in low and high densities. In the previous 
work of Prakash et al [H [9] the parameterization of the potential term of the SE is achieved by the introduction 
of three different choices of the potential contribution to the SE. To advance, in the present work, we suggest a 
more generalized parameterization of the potential term of the SE which is more flexible and efficient, reproduces 
the predictions of most microscopic calculations of dense matter [25| and confirms the results of various empirical 
data. 

The most striking feature of the proposed parameterization is the different form of the parameterization function 
F{u) for densities below saturation point and for densities above this point. This is not surprising since this was 
already entailed in microscopic calculations. Although the behavior of the SE for densities below the saturation 
point still remains unknown, significant progress has been made only most recently in constraining the SE at 
subnormal densities and around the normal density from the isospin diffusion data in heavy- ion collisions [42[ I43j . 
This has led to a significantly more refined constraint on neutron-skin thickness of heavy nuclei [HI US] and the 
mass-radius correlation of neutron stars [33l HU |35] . For densities above the saturation point the trend of the SE 
is model dependent and exhibits completely different behavior. 

The above characteristic of the SE is well refiected in our proposed models. In view of the previous comment, the 
proposed parameterization of the potential term of the SE has the advantage to be able to reproduce microscopic 
calculations in cases where the SE, at low densities, increases along with the density and then begins to fall although 
the density continues to increase. This is a well known characteristic of a class of Skyrme interactions [T^l dZl dH] 
and of Gogny Hartree-Fock calculations [33] |34l [35] . Special effort has been devoted to find analytical relations 
between the radius R and the pressure P which correspond to a special density n for a fixed value of the mass M 
of the neutron star. So an accurate determination of a neutron star radius will permit evaluation of the pressure 
of neutron star matter. All the above will provide a direct determination of the density dependence of the nuclear 
SE at these densities [571liS]. 

Finally, we also address on the problem of neutron star cooling [JT] HH] ■ It is well known that the direct Urea 
process can occur in neutron stars if the proton concentration exceeds some critical value in the range of 11-15 %. 
The proton concentration can be determined by the trend of the SE especially just above the equilibrium density. 
So, the detailed knowledge of the SE behavior is crucial for the existence of the direct Urea process. 

The plan of the paper is as follows. In section 2 the proposed model and the relatives formulas are discussed 
and analyzed. Results are reported and discussed in section 3, while the summary of the work is given in section 4. 



2 The model 



In general, the energy per baryon of neutron-rich matter may be written as 

Eir^ = + (1 - 2xrEilUn) + (1 - 2xfEitl{n) + ■■■ , (1) 

To a good approximation, it is sufficient to retain in the above expansion only the quadratic term. Thus, the above 
equation takes the form 

^ = ^ + (l-2x)^i.il(n), (2) 



where n is the baryon density (n = + rip) and x is the proton fraction (x — Up/n). The symmetry energy 
Esym{n) = Elym{n) can be expre 
and symmetry [x = 1/2) matter 



Esym{n) — Esyrn{n) Can be expressed in terms of the difference of the energy per baryon between neutron {x = 0) 



Esym[n) - — -J ■ (3) 



In the present work we consider a schematic equation for symmetric nuclear matter energy (energy per baryon 
E / A or equivalently the energy density per nuclear density e/n) which is given by the expression [1] 
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where = (hkp)'^ /2mN is the Fermi energy per baryon in equilibrium state and Ug is the saturation density. 

The density dependent potential energy per nucleon V{u) of the symmetric nuclear matter is parameterized, 
based on the previous work of Prakash et. al. [U [9] as follows 
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where pf is the Fermi momentum, related to p^p hy pp = p^u^^^. The parameters Ai and A2 parameterize the 
finite forces between nucleons. The values used here are Ai = 1.5p^ and A2 — Sp'^. The parameters A, B, B' , 
(7, Ci and C2 are determined with the constraints provided by the properties of nuclear matter saturation. In the 
present work the values of the above parameters are determined in order that E{n = ns)/A — rriNC^ = —16 MeV, 
Ug = 0.16 fm^-^ and Kq = 240 MeV. In general the parameter values for three possible values of the compression 
modulus Kq [Kq — 9nQ^-^-^\ng) are displayed in table I, on Ref. [1 . 



To a very good approximation, the nuclear symmetry energy Esym can be parameterized as follows [B] 



Esyra 



{u) = (2^/3 - 1) ^-El _ F{u)) + SoF{u), (6) 



where Sq is the SE at the saturation point, Sq = Esyjn{u — 1). In general, theoretical predictions give 6*0 = 25 — 35 
MeV. In the present work we consider Sq — 30 MeV. The function F{u) parameterizes the potential contribution 
of the nuclear SE and has to satisfy the constraints F{u = 0) = and F{u = 1) = 1. Equation © can be written 
in a more instructive form by separating the kinetic and the potential contribution of the SE. 

E,y,n{u) ~ 13?^+ 17F{u) . (7) 

Kinetic Potential 

In the previous work of Prakash et. al. [1], three representative forms that mimic the results of most microscopic 
models are used and have the following form 

2ii2 

F{u) = u, F{u) = -—, F{u)^V^. (8) 
1 + u 

In the present work we generalize the previous form of the function F(u) in two ways. First, the function F(u) 
is parameterized as follows 

F{u) = (9) 

where the parameter c (hereafter called potential parameter) varies between 0.4 < c < 1.5 in order to get reasonable 
values for the SE. It is obvious that according to the above formula the trend of the potential part is the same 
both in low and high values of the baryon density. 

The information gained from microscopic theoretical calculations shows that this is not the general case for 
the potential part of the SE. On the contrary, the SE exhibits different trends in low and high densities. So, one 
should try to find a new formula for the function F{u) which satisfies the above restrictions. In the spirit of the 
previous statement we propose a new parameterization of the function F(u). The new function which is more 
flexible compared to the previous ones, reproduces the SE for most realistic calculations and has the following form 

u < 1 

F{u) = { (10) 
u'^^ei-" + (u - l)(ci + 1 - C2) u>l. 

The function F{u) satisfies the constraints F{u 1+) = F{u 1") and F'{u 1+) = F'{u 1"). The 
derivative of the function, compared to equation ([5]), is determined by the parameters ci and C2 (hereafter called 
potential parameters). 

In order to construct the nuclear equation of state, the expression of the pressure is needed. In general, the 
pressure, at temperature T = 0, is given by the expression 

„ ^d(eln) de ,,,, 
P = n^ ^ ' = n— ~ e. (11) 
dn dn 



From equations , ^ and (jlip we found that the contribution of the baryon to the total pressure is given by 
the relation 
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(12) 



The leptons (electrons and muons) originating for the condition of the beta stable matter contribute also to the 
total energy and total pressure [6]. To be more precise the electrons and the muons which are the ingredients of 
the neutron star are considered as non-interacting Fermi gases. In that case their contribution to the total energy 
and pressure is given by 



4^5 r 



(2z3 + z)(l + z2)i/2 _sinh-i(z) 



(2z^ - 3z)(l + + 3sinh-^(z) 



(13) 



(14) 



where z — kp/mic. Now the total energy and pressure of charge neutral and chemically equilibrium nuclear matter 
is 

etot = efc + ^ e;, (15) 

l—e~ .j.L~ 
l—e~ ,fi~ 

From equations (fT5|) and ([T6|) we can construct the equation of state in the form e = e(P). What remains is 
the determination of the proton fraction x in /3-stable matter. In that case we have the process 



■ p + e +Ve 
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(17) 



that takes place simultaneously. We assume that neutrinos generated in these reactions have left the system. This 
implies that 

/i = ^„ - = ^e, (18) 

where /x„,/ip and /ig are the chemical potential of the neutron, proton and electron respectivelly. Given the total 
energy density e = e(n„,np), the neutron and proton chemical potential can be defined as 



de 



r\ \il"p 1 

It is easy to show that after some algebra we get 
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dE_ 
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In /3 equilibrium one has 



dE 

9^ = -g^(Eb{n,x)+E,ix)) = 0, 



(19) 

(20) 
(21) 

Up = 

(22) 

(23) 
kp^c 

(relativistic electrons). Equation (|23|) determines the equilibrium proton fraction x{n) once the density dependent 
symmetry energy Esym{n) is known. After straightforward algebra we get 



where Eb{n, x) the energy per baryon and Ee{x) the electron energy. The charge condition implies that rig 
nx or kp^ = kp^. Combining the relations ^ and ([^0]) we get 

/i = 4(1 - 2x)Esym{n). 

Finally by combining equations (fT8|) and ([22|) we arrive at the relation 

4(1 - 2x)Esyrn{n) = ?ic(37r2rie)^/^ = hc{?,'n:'^nxf/'^ , 



where we considered that the chemical potential of the electron is given by the relation /ig = ^/ k'^p c? + m^c' 



2^4 



where 



x[n) = - - - ([2/3(7 - 1)]'/' - [2/5(7 + 1)]'/') , 
/3 = 3^2„(?ic/4ii;,,™(n))3, ^ = A + M") 



(24) 
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When the electrons energy is large enough (i.e. greater than the muon mass), it is energetically favorable for 
the electrons to convert to muons 

— > /i^ + Pf_i + Ve- (25) 

However, in the present work we will not include the muon case to the total equation of state since the muon 
contribution does not alter significantly the gross properties of the neutron stars. 

It is worthwhile to notice that the present model satisfies the relativistic causality. That means the speed of 
sound which was defined from the relation, 

CsV _dP _ dP/dn 
Ci J de de/dn ' 

does not exceed the speed of light for any value of the baryon density. This is a basic treat for any realistic EOS, 
regardless the details of the interactions among matter constituents or the many body approach . 

The most efhcient process, which leads to a fast cooling of a neutron star, is the direct Urea process involving 
nucleons 

n — > p + e~ + Ve, p + — > n + i/e ■ (27) 



This process is only permitted if energy and momentum can be simultaneously conserved [47[ I48j . This requires 
that the proton fraction must be x > 1/9 ~ 0.11. From equation ((24|) it is obvious that the proton fraction x is 
sensitive to the density dependence of the SE and as consequence to the parameterization of the potential part 
of Egym{n). So, in the present work it is worthwhile to study the relation of the proton fraction and the relative 
parameterization and also to check if our parameterization satisfies the constraints for the beggining of the Urea 
process. 

In order to calculate the gross properties of a NS we assume that a NS has a spherically symmetric distribution 
of mass in hydrostatic equilibrium and is extremely cold (T — 0). Effects of rotations and magnetic fields are 
neglected and the equilibrium configurations are obtained by solving the Tolman-Oppenheimer-Volkoff equations 
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dP{r) _ Gm{r)p{r) f P{r) \ / 4:nr^P{r)\ / 2Gm{r) 
dr \ cp-p{r) J \ c^m(r) / \ c^r 

dMir) ^ ^ , ^ inPeir) 

— ^ = Wp(r) = (28) 

dr 

To solve the set of equations (|28|) for P{r) and A/(r) one can integrate outwards from the origin (r — 0) to 
the point r = R where the pressure becomes zero. This point defines R as the coordinate radius of the star. To 
do this, one needs an initial value of the pressure at r = 0, called Pc — P{r ~ 0). The radius R and the total 
mass of the star, AI = M{R), depend on the value of Pc- To be able to perform the integration, one also needs 
to know the energy density e(r) (or the density mass p(r)) in terms of the pressure P{r). This relationship is the 
equation of state for neutron star matter and in the present work has been calculated for various cases by using 
our model. It should also be noted that besides the stellar radius and mass, other global attributes of a neutron 
star are potentially observable, including the moment of inertia and the binding energy |371 149] . Thus, it would be 
of interest to study the nuclear symmetry dependence on these attributes. Such work is in progress. 



3 Results and discussion 



First we apply our model in a simple case where the potential part of the SE is parameterized as F{u) = u'^ and 
the total SE contribution can be written as follows 



(29) 



The potential parameter c varies between 0.4 < c < 1.5 which gives reliable values of the SE. The total pressure 
of the cold beta-stable nucleonic matter is given by 



P(n, x) = 71^ 



A 



+ P' ™(n)(l-2x)^ 



+ Pe-{n), 



where E'^y^in) takes the form 



Egymi'^) 



_ dEsymjn) _ J_ 
dn rt. 
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-1/3 + 17CU--1 



(30) 



(31) 



We are interested for the total pressure at the saturation density rig. Considering that the electron pressure 
Pe- (n) is [37] 

Pe- {n) = nxil - 2x)Esym{n) (32) 
then the total pressure at Ug is given by the expression [37j 

Ps{ns,Xs) = ns{l ~ 2xs) [nsE'gy^{ns){l - 2xs) + Esym{ns)xs\ , (33) 



where E'^y^^iris) 



and the equilibrium proton fraction at is given 
Xs ^ {ZTi''ns)-\AEsyra{ns)/ncf ~ 0.04. (34) 



For small values of we find that 



Ps{ns,Xs)~nlE' (us). (35) 



From the former expression it is obvious that the pressure is mostly sensitive to the density dependence of the 
SE at the saturation point n^. Using our model from equation (j3ip we get 



Ky^{ng)^^(^ + nc]. (36) 



From equations ([35]) and ([36]) it is concluded that the relation between the pressure Pg and the potential 
parameter c is 

P,{n„xs)^n,{^ + n(^ . (37) 

In order to calculate the global properties of the neutron star, radius and mass we solved numerically the 
TOV equations (|28p with the given equations of state constructed with the present model. For very low densities 
{n < 0.08 fm"'^) we used the equation of state taken from Feynman, Metropolis and Teller [5^ and also from 
Baym, Bethe and Sutherland [51]. 

In order to illustrate the density dependence trend of the SE proposed in our model we display in figure la 
Esym{n) as a function of the density n for various values of the potential parameter c. It is obvious that the 
parameter c affects decisively the trend of the SE, especially at high values of the density. So, it is very interesting 
to study how the values of the parameter c, and consequently the potential contribution, affect the gross properties 
of the NS. In the same figure, results of Ref. [39] (the case Al8+(5u+UIX*, see TABLE VI and VII of Ref. 39J) 
are included. It is found that the use of the phenomenological equation ([7|) with proper value (c = 0.9) of function 
^ reproduces the results of the above microscopic calculations. 

Figure 2a demonstrates the linear dependence between the radius Rmax and the parameter c. From our analysis 
it is concluded that there is a direct relation between Rmax and the parameterization of the potential part of the 
SE. The star symbol corresponds to the case A18+^u+UIX* with ci — 0.9. The corresponding relation was derived 
with the least-squares fit method and has the form 

R,nax = 9.10195 + 2.08304c. (38) 

In addition, in figure 2b we indicate the behavior of the neutron star radius Ria, which corresponds to a 
neutron star mass M — 1.4Mo, versus the parameter c. The star symbol corresponds to the case A18-|-(5u-|-UIX* 
with ci = 0.9. It is obvious that there is also a linear relation between i?i.4 and c which is 

RiA ^ 10.52114 + 3.56746c. (39) 

Figure 2c displays the correlation between the maximum mass of the neutron star M^ax and the parameter c. 
The star symbol corresponds to the case AlS+Ju+UIX* with ci — 0.9. We found an almost linear relation between 
Mmax and c which has the form 

Mmax = 1.87442 + 0.12344c. (40) 

By combining equations ([36| and (|39|) we found a linear relation between Ria and E'^y^ijis) which has the 
form 

RiA = 8.702 + amZQE'gy^in,), (41) 

and vice- versa 

KymM = -259.185 + 29.783i?i.4. (42) 



n 



In order to illustrate further the relation between the radius Ria and the trend of the SE we plot in figure 3a 
the radius Ria versus the derivative of the symmetry energy E'gy^{2ins/2) at the baryon density n — 3ns/2. A 
linear relation is found, which has the form 

RiA = 10.1798 + 0mA2E'^y^{3ns/2). (43) 

It is concluded that there is a direct relation between the radius Ria and the trend of the SE, close to the 
saturation point rig. In addition in figure 3b we plot E'^y^{3ns/2) versus Ria with the linear correlation 

Kym{^n,/2) = -420.40895 + 41.3066i?i.4. (44) 

To ensure the relativistic causality in the present model we display figure 4a, where the ratio of the speed of 
sound to speed of light Cg/ci is plotted versus the baryon density for various values of the parameter c. Evidently, 
in all cases the speed of sound does not exceed that of light even at high values of the baryon density. In addition 
in figure 4b we plot the proton fraction Xp calculated from expression ([M]) as a function of the baryon density n. 
It is obvious that only in the cases where c > 0.5 the proton fraction, after a specific density, exceed the critical 
value x^^'^'^ ~ 0.11 which ensures the beginning of the Urea process. 

We also tried to find the correlation between the pressure P (and consequently the radius R) and the SE for 
other values of the density n. In order to clarify the problem of the expected relation between the radius and the 
pressure we present a more simplified model of a non-relativistic equation with a polytrope type of EOS. Thus the 
EOS has the form [3 [37] 

P^Kp^, ^=1 + 1. (45) 
A 



and the radius of the star is given by 



R = 



{\ + l)K 



(46) 



where pc is the central density and is the solution of the equation 0(^i) — 0, where the function 6'(^) is the 
solution of the differential Lane-Emden equation 

Now it is obvious from equations (|45|) and ([46]) that in case A = 1 (or 7 = 2) where ^ = tt we have 



R 

pl/2 



2G 



1/2 ^ 



(48) 
P 



Thus, from equation ()48|) we concluded that in the case of a polytrope with 7 = 2 there is a universal relation 
of the ratio RjP^I'^ calculated for a specific value of the density p. However if general relativity effects are included 
in the above analysis the exponent 1/2 of the pressure is found to be smaller [37'. 

Following the above statement we plot in figure 5a the quantity RiaP^"" as a function of the radius R\a- One 
can see that there is a correlation between the radii R\a and the pressure evaluated at densities Iris, ins/ 2 and 2ns. 
The values of the parameters a and C(n) have been defined by least-squares fit of the expression Ria ~ C{n)P°' . 

In addition, in figure 5b we plot the quantity RiaR^^^^ as a function of the Ria to compare it with the previous 
work of Lattimer et. al [37| . It is worthwhile to notice that the quantity RiaR^^^'^ is a mild increasing function 
of the radius Ria- This effect is more evident for densities far from the saturation( n = 3ns/2, 2ns). So, from our 
study it is concluded that there is a slight dependence of the quantity RiaR^'^^'^ from the potential parameter c 
and consequently from the trend of the SE. 

We proceed now in the more complicated case where the function F{u) is given by expression (llOp . In that case 
the derivative of the Egym (n) is given by 



j_ r26„ 

Us L 3 

Kymin) = { (49) 
J_ [f u-1/3 + 17 + 1 _ C2 + gi-^w'^^ (f - 1))] U>1 . 

The potential parameters ci and C2 varied between 0.5 < ci < 1.2 and < C2 < 2 in order to get a reliable 
density dependent SE. 

In figure lb we display Esym (n) as a function of the density n for various values of the potential parameters ci 
and C2. In general the case is as follows, for fixed values of the parameter C2, the SE is an increasing function of ci. 



In addition, for fixed values of the parameter ci the increase of the parameter C2 leads to a decrease of the SE. It is 
seen that, within the present model, the stiff or soft behavior of Esymin) found in various microscopic calculations, 
is reproduced. As a comparison, similar to figure la, results of Ref. [39] (the case A18+(Su+UIX*) are included. It 
is found that the use of the phenomenological equation ([7]) with proper values (ci = 0.77 and C2 = 1.09) of function 
(fTO|l reproduces the results of the above microscopic calculations. 

Figure 6a illustrates the behavior of the radius Ria as a function of the second potential parameter C2 for 
various values of the first potential parameter ci . The calculated points for various values of ci can be reproduced 
by a second order polynomial. 



RiA = 12.58956 + 0.05378c2 - 0.46172c^, ci = 0.5 

RiA = 13.04786 - 0.02752C2 - 0.32340c^, ci = 0.7 

RiA = 13.85705 - 0.08087c2 - 0.21393c^, ci = 1.0 

RiA = 14.61946 -0.19441c2 - 0.14536ci ci = 1.2 . (50) 

In all examined cases, the radius Ria is a decreasing function of the potential parameter C2. This is a direct 
consequence of the softening of the equation of state due to increase of the parameter C2. 

In addition, in figure 6b the behavior of the radius Ria as a function of the first potential parameter ci is 
reproduced for various values of the second potential parameter C2 . The least-squares fit values are given for the 
following linear equations 

RiA = 11.09603 + 2.85172ci, C2 = 0.0 
RiA = 11.01810 + 2.85517ci, C2 = 0.5 

RiA = 10.42034 +3.06724ci, C2 = 1.2 . (51) 

Unlike the previous case, Ria is an increasing function of the potential parameter ci. The increase of the 
parameter ci leads to the stiffness of the SE as indicated in figure 6b. It is worthwhile to note that the slopes of 
the best fit lines are almost the same and there is just a shift of the lines depending on the values of the parameter 

C2- 

Also, from figure 6a and 6b we conclude that the radius Ria depends mainly on the parameter ci which 
determines the derivative of the Esym{n) and also the pressure Pgat at the saturation density Ug. However, there 
is a small dependence on the parameter C2 which is connected with the trend of Esymin) at higher values of the 
density Hs . Figures 6c and 6d demonstrate the dependence of the radius Rmax and the mass Mmax respectively on 
the parameter ci (for fixed values of the parameter C2). The star symbol corresponds the the case A18+(5u+UIX* 
(ci = 0.77 and C2 = 1.09). In both cases a linear relations holds between Rmax, Mmax and the parameter ci. The 
lines correspond to the least-squares fit values. 

It is of interest to compare the NS properties {Ria, Rmax and Mmax) which originated from the use of equations 
^ and (fTO|) . As an example we use the parameterization of equations ([9]) and (jlOp which reproduce very well the 
Esym trend of the case A18-|-(5u-|-UIX*. As a result it is found that Ria{c = 0.9) = 13.42 and Ria{ci = 0.77, C2 = 
1.09) = 12.76 (difference 5 %) , Rmax{c = 0.9) = 11 and Rmax{ci = 0.77, C2 = 1.09) = 10.59 (difference 3.7 
%), Mmaxic = 0.9) = 1.978 and Mmaxici = 0.77, C2 = 1.09) = 1.974 (difference 0.2 %). It is obvious that in 
the case Al8-f (5u-fUIX* we receive almost identical results for Mmax while there is a small difference for Ria and 
Rmax- The differentiation of the values of Ria is a consequence of the linear relation which hold between the radius 
RiA and the derivative of the Esym (and consequently according to ([55)) to the pressure), close to the saturation 
point (see figures 3 and 7). More specifically, we receive for the two cases, Ps{c = 0.9) = 3.83467 MeV fm~'^ and 
Psici — 0.77, C2 = 1.09) = 3.48107 MeV fm^'^. In general, the small differentiation on the radii is not surprising 
since the trend of equations ([9]) and (|10p. due to suitable parametrization, are similar. However, it is worth to point 
out that equation pOjl is a generalization of equation ([9]) , in the meaning that while equation ([9]) describes well the 
case where Egym is a increasing function of the density, equation ([TOj) is sufficiently flexible to describe in addition 
the case where Esym at low densities increases along with the density and then begins to decreases although the 
density continues to increases. 

To illustrate further this point, we studied the correlations between the derivative of the symmetry energy 
E'sym and the radius Ria close to the saturation point n = 'ins/2. In figure 7a we plot the radius Ria versus the 
derivative of the symmetry energy E'symi^T^s/^.) for fixed values of the potential parameter C2. One can see that 
there is a linear relation between Ria and E'^y^{3ns/2) just like in figure 3a. The effect of the parameter C2 is to 
induce a parallel shift of the best fit Hues. In figure 7b we indicate the inverse relation, that means E'^yj^{3ns/2) 
versus Ria- The least-squares fit values for both cases and for various values of the parameter C2 are given for the 
following equations 

RiA 8.70394 + 0.02684£^j,„(37is/2), C2 = 



E'sym{^ns/2) = 



-318.5026 + 36.82989i?i.4 , 



(52) 



RiA = 9.73291 + 0.02687£;^j,„(3ns/2), 
Kym{^nj2) = -359.07613 + 36.98138iii.4 , 



C2 = 0.5 



(53) 



RiA = 10.27305 + 0.02887£;^y„(3n^/2) 
E' {3nj2) = -354.58213 + 34.54212i?i.4 . 



= 1.2 



(54) 



In figure 8a, likewise with the figure 4a we display the ratio c^/c; as a function of the density for various cases. 
It is obvious that the relativistic causality is satisfied once again. In figure 8b we display the proton fraction Xp 
as a function of the density for various cases. A more systematic study of the Xp leads to the conclusion that the 
potential parameter ci plays the most critical role for the occurence of the Urea process. Specifically a higher value 
of the ci leads to the beggining Urea process in smaller values of the baryon density. 

Figure 9 illustrates the behavior of the quantity R\aP~°' as a function of the radii i?i.4 for pressure determined 
at n = n^, 'ins/2, 2ns, and also for C2 = (figure 9a), C2 = 0.5 (figure 9b) and C2 = 1.2 (figure 9c). 

In figure 10 we plot the quantity RiaP'^^"^ as a function of i?i.4 for the pressure determined at n = rig, 3/1^/2, 
2ns and for C2 = (figure 10a), C2 = 0.5 (figure 10b) and C2 = 1.2 (figure 10c). It is obvious once again that 
the quantity RiaP~^^^ is almost constant only when the pressure is calculated at the saturation point n.,. When 
the pressure is calculated at densities n = 3ns/2 and n = 2ns the quantity RiaP~^^^ is an increasing function 
of the radius R\a- Thus, as in the case of the simple parameterization of the SE, it is concluded that there is a 
dependence of the quantity RiaP~^^^ from the first potential parameter ci as well as from the second potential 
parameter C2 and consequently from the trend of the SE both for low and high values of the baryon density. 

4 Summary 

In the present work we performed a systematic study of the effect of the potential part of the SE on the global 
properties of neutron stars (masses, radii and composition). The potential part of the SE was parameterized in 
a generalized form both for low and high values of the baryon density in order to be efficient in reproducing the 
results of most microscopic calculations of dense nuclear matter. 

In the case of the simple parameterization of the SE the most striking feature of our study was the derivation 
of a linear relation which stands between the maximum mass M^ax, the radius Rmax and the radius Ria with 
the potential parameter c. In addition, a linear relation stands between the i?i.4 and the derivative of Esy^{n) for 
densities close to the saturation point (n = ns,3ns/2). It was concluded that quantity R\aP~°' (with a and C{n) 
fitting parameters) appears to be constant for the densities n = n^, 3ns/2, 2ns. However, the quantity RiaP~^^'^ 
exhibits an increasing behavior as a function of Ria for n — 3ns/2, n — 2ns. 

In the case of the more complicated parameterization, where the SE is parameterized in a different way for low 
and high values of the density, similar results were taken. Specifically i?i.4 is a function of both potential parameters 
ci and C2. This means that the value of i?i.4 is affected from the density dependent trend of the SE, both in low 
and high densities. However, we showed that for fixed values of the parameter C2, close to the saturation point 
(n = 3ns/2), a linear relation between the Ria and the E'sy^{3ns/2) stands. Finally, the quantity RiaP~'^ = C{n) 
appears to be constant after a suitable parameterization of the parameters a and C (n) but still remains dependent 
from the second potential parameter C2. The quantity RiaP~^^'^, as in the previous case, exhibits an increasing 
behavior as a function of the R\a for density values above the saturation point. 
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Figure 1: (a) Esym{n) for varfous values of the potential parameter c of the function F{u), given by the equation 
((HI), versus the baryon density n. The case A18+(5u+UIX* from Ref. ,39j is also included with the least-squares fit 
curve which corresponds to the value c = 0.9 (b) Esymiji) for various values of the potential parameters ci and C2 
of the function F{u), given by the equation (fTO|) . versus the baryon density n. The case A18-t-(5u+UIX* from Ref. 
[55] is also included with the least-squares fit curve which corresponds to the values ci 0.77 and C2 = 1.09. 




Figure 2: (a) The radius Rmax as a function of the parameter c. (b) The radius Ria as a function of the parameter 
c. (c) The maximum mass Mmax of the neutron star as a function of the parameter c. The solid lines correspond to 
the least-squares fit expressions (LSF) (l38t . ((39)) and (I40p respectively. In all figures the star symbol corresponds 
to the case A18 + 5u + UIX* ( c = 0.9). 
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Figure 3: (a) The radius Ria versus the derivative of the Symmetry Energy i?^^^(3n<,/2). (b) i?^y„(3ns/2) versus 
RiA- The sohd Hues correspond to the least-squares fit (LSF) expressions (|43)) and (|44|) respectively. 




Figure 4: (a) The ratio Cg/ci versus the baryon density n for various values of the potential parameter c. The 
line for c = 0.9 corresponds to the case A18 + 6u + UIX*. (b) The proton fraction Xp versus the density n for 
various values of the potential parameter c. The line for c — 0.9 corresponds to the case A18 + Su + UIX*. The 
short-dotted line shows the beginning the direct Urea process (xp = 0.11). 
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Figure 5: (a) The quantity i?P^° as a function of the radius Ria for pressure determined at n = n^, n — 3n 
and n = 2ns ■ For each density, the least-squares fit value for the exponent a is indicated, (b) The quantity RP~ 
as a function of the radius Ria for pressure determined at n = Ug , n = 3ns/2 and n = 2ns. 





Figure 6: (a) The radius Rn as a function of the second potential parameter C2 for various values of the first 
potential parameter ci. The lines correspond to the least-squares fit expressions ()50p . (b) The radius Ria as 
a function of the first potential parameter ci for various values of the second potential parameter C2. The lines 
correspond to the least-squares fit expressions ([51]) . (c) The radius Rmax as a function of the first potential 
parameter ci for various values of the second potential parameter ci. The lines correspond to the least-squares 
fit expressions, (d) The maximum mass M„iax for various values of the second potential parameter ci. The lines 
correspond to the least-squares fit expressions. In all figures the star symbol corresponds to the case Al%+5u+UIX* 
( ci = 0.77 and C2 = 1.09) 




Figure 7: (a) The radius Ria versus the derivative of the symmetry energy E'gy,^{3ns/2) for various values of the 
potential parameter C2. (b) E'^y^(3ns/2) versus Ria- The lines correspond to the least-squares fit expressions 
((52| . ((53| and ^ respectively. 
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Figure 8: (a) The ratio c^/q versus the baryon density n for various values of the potential parameters ci and 
C2- The line for ci = 0.77 and C2 = 1.09 corresponds to the case A18 + 6u + UIX*. (b) The proton fraction Xp 
versus the density n for various values of the potential parameters ci and C2. The line for ci = 0.77 and C2 = 1.09 
corresponds to the case A18 + 6u + UIX*. The short-dotted line shows the beginning of the direct Urea process 
(xp = 0.11). 
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Figure 9: (a) The quantity i?,_P^° as a function of the radius Ria for pressure determined at n = Ug, n = ius/i 
and n = 2ns for C2 = 0. (b) The same as before for C2 — 0.5. (c) The same as before for C2 = 1.2. For each density, 
the least-squares fit value for the exponent a is indicated. 
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Figure 10: (a) The quantity RP as a function of the radius R\a for pressure determined at n = rig, n = 'ins/2 
and n = 2ns for C2 = 0. (b) The same as before for C2 = 0.5. (c) The same as before for C2 = 1.2. 



